Vortices in magnetically coupled superconducting layered systems 
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Pancake vortices in stacks of thin superconducting films or layers are considered. It is stressed 
that in the absence of Josephson coupling topological restrictions upon possible configurations of 
vortices are removed and various examples of structures forbidden in bulk superconductors are given. 
In particular, it is shown that vortices may skip surface layers in samples of less than a certain size 
Re which might be macroscopic. The Josephson coupling suppresses Re estimates. 
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I. INTRODUCTION 



Attempts to measure the anisotropy parameter in Bl- 
and Tl-based high- Tc superconducting compounds yield 
very large values ~ 100 — 300.EI In some layerpjd organic 
superconductors this parameter is even highcrfl This has 
led many to believe that the Josephson coupling in these 
materials is so weak that for many purposes it can be 
disregarded altogether. Models of vortices in these com- 
pounds based on pure magnetic coupling between two- 
dimensional pancake vortices proved to suffice for many 
applications .ffEl 

We consider in this paper a system of thin supercon- 
ducting layers coupled only via the magnetic field be- 
tween them. This system is qualitatively different from 
the bulk superconductors: there is no phase coherence 
across the layers. In three-dimensional bulk materials 
vortices are banned from terminating inside supercon- 
ductor because of topology: the phase changes by 27r 
when one circles the vortex core at which the phase is 
singular. A core termination would have meant that the 
phase could acquire the 2-n change along a path that does 
not contain a singularity. 

Clearly, this ban does not hold in layered systems with 
no Josephson coupling. A vortex perpendicular to the 
layers or, better to say, a stack of pancake vortices (or 
"pancakes"), may terminate, "'^ in principle, at any layer 
and channel the flux at least partially into interlayer 
space. This event might be energetically unlikely, but 
it is not forbidden by topology of the phase, which is 
defined separately in each layer. 

We argue below that within the model of layered ma- 



terials with zero Josephson coupling such termination in 
layers adjacent to the surface might be energetically fa- 
vorable in samples of finite size. We use the formal tech- 
nique suggested by one of us,B which is reviewed briefly. 
The technique allows us to compare energies of various 
configurations in a straightforward manner. As one of the 
applications, we compare energies of two configurations: 
one with a vortex piercing all layers of a half-space mul- 
tilayer and another with the vortex terminating in one 
of the top layers. We show that the energy cost of the 
subsurface vortex termination diverges with the sample 
size R: 
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where A = A^^/s, \ab is the penetration depth, s is the 
layer periodicity, k = Xab/^ab-, and ^ab is the coherence 
length. The divergence of Ae is weak, however, and in 
samples smaller than Re '^ nXab the subsurface termina- 
tion becomes energetically favorable. With Xab ~ 0.3 /im, 
K « 100 for Bi-2212, we estimate Re ^ 30 fim with a large 
numerical factor so that Re may reach a macroscopic size. 
This estimate is reduced considerably by including the 
Josephson coupling because one has to include in the bal- 
ance the energy of Josephson strings channelling the flux 
sideways into the interlayer space. We estimate, however, 
that even then Re remains on the order of microns. This 
suggests that if the surface of a layered compound like 
Bi-2212 has defects separated by Re, the occurrence of 
vortex cores at the surface is a rare event. We also show 
that for a vortex terminating at a depth <C Xab under 
the surface, most of the magnetic flux crosses the core- 
free layers into the free space, the absence of the core 



'^We use the term "vortex termination in the layer N" for 
the situation when the vortex cores are present in the layer 
N and in all layers under, whereas the cores are absent in the 
layers above N up to the the sample surface. 



notwithstanding. Hence, we expect vortices terminating 
under the surface to be invisible for the scanning tunnel- 
ing technique, but detectable in decoration experiments. 
Other examples of possible but unusual configurations 
of pancake vortices in layered materials are considered. 
It is shown that a vortex in a film of a finite size and 
containing several superconducting layers may have nor- 
mal cores only in internal layers and carry flux different 
from the flux quantum (J)q . The last feature appears also 
in exotic configurations such as a stack of pancake vor- 
tices situated in every other layer. Implications of these 
possibilities are discussed. 



II. APPROACH 

We begin with a brief reviewi-ef thin superconducting 
films. As was stressed by Pearl,t3 the situation in a thin 
film differs from that of a bulk since a large contribution 
to the energy of a vortex comes from the stray fields. In 
fact, the problem of a vortex in a thin film is reduced 
to that of the field distribution in free space subject to 
certain boundary conditions at the film surface. Since 
curlh = divh = outside, one can introduce a scalar 
potential for the outside field: 



V(p, vv = o. 



(2) 



To formulate the boundary conditions for the outside 
Laplace problem let us consider a film of thickness d ^ X 
occupying the xy plane; A is the bulk penetration depth 
of the film material. For a vortex at r = 0, the London 
equations for the film interior read: 
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(3) 



where z is the unit vector along the vortex axis. Averag- 
ing over the thickness d we obtain 
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(4) 



where g(r) is the sheet current density, r = {x,y), and 
A = X^ /d is the film penetration depth. Other compo- 
nents of London equation turn identities after averaging. 
Since all derivatives d/dz are large relative to the tan- 
gential d/dr, the Maxwell equation curlh ~ 47rj/c is re- 
duced to conditions relating the sheet current to discon- 
tinuities of the tangential field: 



in 



-9x 



K-- 



An 



-9y 



ht 



(5) 



Here the superscripts ± stand for the upper and lower 
[z — ±(i/2) faces of the film. The field component per- 
pendicular to the film, /i^, is the same at both film faces. 
Substituting Eq. (^) into Eq. (^) and using divh = 0, 
we obtain: 
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dz dz 

This equation along with 



= 0o^(r) . 



(6) 



(7) 



and conditions at infinity constitute the boundary con- 
ditions for the Laplace problem, Eq. (0), of the field 
distribution outside the film. 

Let us turn to the question of energy. We consider a 
general situation of vortices in finite bulk samples. The 
energy consists of the London energy (magnetic -I- ki- 
netic) inside the sample, e'*', and the magnetic energy 
outside, e^"': 
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Then, for the potential gauged to zero at cx) (which is 
possible in zero applied field) one has 
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where the integral is over the sample's surface with dS 
directed inward the material. If (as it may happen in 
some situations considered below) (p does not vanish at 
infinity, one has to examine the integral (^ at infinity. 

The London part can be transformed integrating the 
kinetic term by parts: 
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where the integral is over the samples surface and the 
surface of the vortex core; there might be more than one 
sample in the system, while not all of them may contain 
vortices. The integral over the samples surface is further 
transformed: 

* dS • (j X Vip) ^ (h dS- ip{V X j) 

(take a closed contour at the sample surface, consider the 
total sample surface as made of two pieces supported by 
this contour, and apply the Stokes theorem to the inte- 
gration over each piece to show that ^ dS • (V x ipj) = 0). 
Combining the result with e'"^) of (H), one obtains 
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The expression in parentheses is zero for samples with 
no vortices, whereas for those containing vortices it is 
(/<oV(5^^-' (r — Ty) where v is the direction of the vortex 
crossing the surface at the point r„ ((5'^'(r — r„) is the 
two-dimensional 6 function). We then obtain 
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^<f dS.(hxj), (11) 

90C J core 



with Tent and r^x being the positions of the vortex entry 
and exit at the sample surface (the vortex is assumed 
to cross the sample surface at right angles). Note that 
if there are more than one superconductor present, but 
the vortex pierces only one of them, the result ( pi] ) still 
holds (although (p's differ for each particular configura- 
tion). For more than one vortex in the system, one has 
to sum up over all vortices. For thin films, the integral 
over the core surface (oc d) can be neglected in Eq. (11). 
It is— jiiistructive to see now how the known Pearl 
resultaljO for vortex in a film and a bulk half-space can 
be obtained within the approach outlined here. We also 
demonstrate an added advantage of the method, a rela- 
tively simple way to evaluate energies. 



A. Vortex in a thin film 

Consider a thin film situated at z = 0. The general 
form of the potential which vanishes at z ^ +00 of the 
empty upper half-space is 



(p(r,z>0) = 
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with k 



/kl + k'^. In the lower half-space we have to 

replace z — > —z in Eq. (p^. The two dimensional (2D) 
Fourier transforms (^„(k)and (pi(k) for the upper and 
lower half-spaces are obtained with the help of boundary 
conditions (6|J3): 

- kipu + Ak'^{(pi - ifu) = 4>o , -kipu = kipi . (13) 

This system yields: 
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A:(H-2fcA) ■ 
We point out first that the total flux crossing a plane 



z = const IS 

$^ = / h^,(Pr = hJk 
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for any z, i.e., the film is crossed by the flux c/jq. And 
the second: according to (111]), the energy e of the Pearl 
vortex is given by 
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where the cutoff at kmax ~ Stt/^ is introduced to a loga- 
rithmically divergent integral. This yields 
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B. Vortex in a half-space 

Let now the half-space z < be filled with a super- 
conductor having the penetration depth A. The stray 
field in the free space z > is given by the potential 
(|l^). Within the superconductor we have the London 
equation (m. The general solution can be written as 



h-h^o^+ht"), 



(18) 



where h'^'^) solves Eq. (g) with zero right-hand side 
(RHS), whereas h^") is a particular solution of the full Eq. 
(^). The latter can be taken as the field of an infinitely 
long unperturbed vortex along z; this assures correct sin- 
gular behavior at the vortex axis. The Fourier transform 
of this field is 



h('^) = 
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(19) 



We now Fourier transform the homogeneous Eq. (ra) for 
h] (r) with respect to the variable r and obtain 

hf'\k, z)(l + e\^) - \^dlhf\\i, z) = , (20) 

where i = x,y, z. The solution which vanishes deep in 
the superconductor is 



hf'\k,z)^H,{k)e'> 
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= VT+A2fc2/A. (21) 



Here, Hi{k)'s are not independent: divh*^°) — gives 



ikxHx + ikyHy + qH^ ^0 . 
The boundary conditions at z = read: 



(22) 



-A:v,„ = i/, + 4"). 



(23) 
(24) 

The four Eqs. ( p2[ - |2^ ) suffice to determine the four un- 
knowns, (p,i and Hx^y.z- We obtain 
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We do not write down Hx,y which describe the spread- 
ing of the vortex field under tbe surface. Equations ( [25| ) 
coincide with Pearl's solution.Ell 

To evaluate the energy of the vortex in this case with 
the help of the general result ( pT| ) , we first calculate the 
potential at the vortex exit: 
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The energy associated with the core, i.e., the integral 
over the core surface in Eq. (O), is 
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Integrating the London relation 
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over a circle of a radius r, we obtain for the current den- 
sity near the core 
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which results in 






(29) 



(30) 



Using Hz of Eq. (|25D with Eqs. (|ll|) and (g6D we obtain 
the energy cost of vortex termination at the surface as 
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where e/, is the energy per unit length of an unperturbed 
vortex. 



III. SUBSURFACE TERMINATION 

A. Half-space + a thin layer 

We now turn to a stack of pancake vortices perpen- 
dicular to the layers of a half-space layered sample. The 
layered structure has the period s whereas each layer is 
characterized by the Pearl length A = A^/d with d being 
the layer thickness. Note that the measurable length Xat 
is related to other lengths in the problem by 
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Let us consider first the stack with a missing pancake 
in the top layer. We model the rest of the stack as a cort. 
tinuous half-space z < with the penetration depth XabSi 
A thin film with the Pearl length A is situated at z = s. 
In the two domains of free space, z > s and < z < s, 
the 2D Fourier transforms of the potential are: 



^«(k)e-'=^ 


z > s , 


(33) 


(pl(k)e'=^+(/.2(k)e-'=^ 


< z < s. 


(34) 



On the film at z = s, we have the boundary conditions 
W) and m). Three more equations are provided by the 
field continuity at z = 0. The requirement divh^''^ = 0, 
Eq. (P4), completes the system thus providing enough 
equations for ipu, ^1,2 and Hx^y^z- 
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D = q-k-{l + 2kA){k + q)e 
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where q is defined in Eq. (|2l|); H^^y is not needed here. 
The fluxes through planes z = const are: 
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z < 0. 



(37) 
(38) 



Within the layer < z < s, the flux $^ decreases from 
the bulk value 00 to the value (p7|). Therefore, a small 
fraction of the flux, 



4>f)Xab/A = 4>QS/Xab ; 



(39) 



is channelled aside between the half-space and the top 
film. 

Using the same formal scheme, we can solve the prob- 
lem when the film at z = s contains a pancake vortex. 
The only difference with the case of "half-space -I- an 
empty film" is in the London boundary condition (M) on 
the film at z = s. We skip details and provide the result 
denoting the new solutions with a star: 
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The energy cost of vortex termination, i.e., the differ- 
ence between the situations without and with the top 
pancake, is given by 



^Ae = [-95e.(0)] - iv^Lti^) - ^L(s) - ^L(O)] 
+ / dz{hz -K): 



(41) 



the last line is the difference in core contributions of vor- 
tex lines in the low half-space, see Eq. ( pOl) 

Evaluating the Fourier integrals for ip's, one notes 
that the integration over k is done within the domain 
2tt/R , 27r/5 where R is the sample size and ^ is the co- 
herence length (usually called ^ab)- For s < £, one can 
expand e~2'=« in the denominator D and see that the 
term with s can be neglected. Indeed, by introducing a 
new variable sinhw = kXab and utilizing A 3> Xab, we 
obtain 



D 



-2k[l + iA/Xab)e' 



-2ke"A/Xab- 



(42) 



This results in a factor 1/A = s/A^^ in front of the inte- 
grals. Then, in the linear approximation in s, we can set 
s = in the integrand after which all integrals become 
simple. We obtain after a straightforward algebra: 
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It is of interest to obtain the current distribution in 
the top film. To this end, we write Eq. (y_2|) in cylindric 
coordinates: 
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The tangential fields above and under the film are: 
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The sheet current is given by 
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where the boundary condition ([7|) at z = s, ipu + fis." 
1^2 = 0, had been used. With the hel p o f Eq. ( p2| ) for D, 
we obtain after integration (see Ref. |I^, 6.663): 

where Ii,Ko are Modified Bessel functions. Thus, 5^ 
vanishes as rla(2Xab/r) when r — » 0, and as I/r for 

r > Xab- 

We note that the same formal scheme can be applied 
to consider the termination of the pancake stack in the 
second, third, etc, layer from the top. The length A 
characterizing the top layers with no cores should then 
be taken as A/2, A/3, etc. 



B. Vortex in a finite stack of layers 

As an example of such a system we consider a pan- 
cake in middle layer of a "short stack" of three films. 
Let the pancake sit at the origin of the film at z = 0; 
the film is characterized by the length A. At z ± s two 
other films are situated having the film penetration depth 
Ai ^ A. Following the same method we write the 2D 
Fourier transforms: 
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(Pci e''^ -I- (pc2 e ''^ , -s < z < . 



Due to the symmetry, hbz{z) = hcz{—z), and Lp^i — 
—ipb2, fc2 — —^bi] this leaves only three coefficients gs 
to be determined. The boundary conditions (|j) and (g) 
at z = s give 
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(49) 
(50) 



The third equation is provided by the London boundary 
condition (ra) at z = 0: 

^bi{l~2kA)~'fb2{l + 2kA) = (t>o/k. (51) 

The continuity condition M) for ft,^ at z = is satisfied 
identically. Equations (^51) yield 
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The fluxes are 

$^(z > s) =0, 
$2(0 < z < s) ^ 4) 
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(l + 2fcA)(l + 2fcAi). 
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Note that the flux (/)os/(2A + Ai) <C 4>o is deflected into 
the space between the layers. For A ;» Ai ::^ s. 
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Note also that for Ai = A 3> s, the flux piercing the 
sandwich made of identical films is 4>o/?>. The energy is 
given by 
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Using the same argument as that leading to Eq. (|4 
we obtain for the denominator 
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Integration in Eq. (|55| ) is now easy. 

Consider the case Ai = A and set s = (everywhere 
except in A = A^/s). Integration yields for the energy of 
one pancake in a stack of three films 
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where it was assumed that the sample size R 3> 27rA and 
the subscript 010 is to indicate that the pancake is situ- 
ated in the middle film. This is to be compared with a 
vortex piercing all three layers, which can be considered 



as a Pearl vortex in a single film with the effective film 
penetration depth A/3. Then, according to Eq. ([l7|), 
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and the subscript 111 indicates that each layer has a pan- 
cake. The difference of energies (^^ and (^8|) is negative 
if 



R < A(A/e)^ 



(59) 



and we have omitted a large numerical factor. This sug- 
gests that in the stack of three films, vortices carrying 
the flux 00 / 3 with a core only in the middle film might 
be energetically preferable to standard vortices piercing 
the whole stack with the flux 0o- 

It is worth noting that because our solution holds for 
A ^ Ai, in fact, it applies also to a vortex piercing N 
middle layers of the 2Ni + N stack of layers provided 
that both Nis and Ns are less than Xab', this allows us 
to treat the Nis as one film with Ai = X^i^/Nis. 



IV. INFINITE DILUTE STACK OF PANCAKES 

We consider now an infinite stack of layers. Let the 
layer at z = and all even layers be free of pancakes; 
the pancakes are situated at z = ±s and at all odd films 
(at r = 0). Consider three regions: (a) s < z < 2s, (b) 
< z < s, and (c) — s < z < 0. 



Z 

2s 
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-s 



4^ 



FIG. 1. The unit cell of a periodic structure < z < 2s. 



The 2D Fourier transforms: 

fa = 'Pal e^" + fa2 6^''^ 
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Due to the symmetry with respect to z = 0, hbz{z) = 
hcz{-z). This gives (pci = ~fb2, fc2 = -fbi- Further 
relations are provided by periodicity which can be ex- 
pressed as haz{z) = hcz(,z — 2s) and yields 



Val 



-'fb2e 



-2ks 



'P>a2 



-Pbie 



2ks 



(63) 



Thus, out of 6 unknown coefficients of the system (pOM62) 
we are left with only two. As the two needed equation 



one can take the London conditions 

z = s: 
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The flux at any plane z =const within the domain (b) 
is readily evaluated: 



^b{z) = hbzik = 0) = [kipbi e""" - kipb2 e ''%=a 
, 2A + Z 6n 
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Similarly, in the domain (c): 
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Therefore, $ k, 0o/2 everywhere throughout the system. 
Since we have in this case only one pancake per the 
period 2s of the structure, the line energy of the stack 
e = ep/2s where tp is the energy per pancake: 
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As above, within the integration domain we can sim- 
plify the denominator, D2 ~ 8fcA, and set s = in 
the integrand. Then we obtain after integration within 

27r/i?, 27r/^: 
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Comparing this with the line energy of a standard vor- 
tex (0o/47rAafc)2 InK we see that dilute stacks considered 
here are possible in small samples (whiskers) of a size 

R < XabK^- 



V. DISCUSSION 

We now estimate how the energy cost (^) of the sub- 
surface vortex termination is affected by the Josephson 
coupling. Clearly, this coupling breaks the cylindrical 
symmetry of the field associated with the straight stack 
of pancake vortices, and the fiux (p9) will not be spread 
even in all directions in the space between the top layer 



and the rest. Still, for small samples R < Xj, in which the 
string cannot fully develop, the asymmetry can be dis- 
regarded. Since the phase difference is zero for perfectly 
aligned stack of pancakes for all pairs of layers except the 
top one, the Joscphson energy can be estimated as 
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On the other hand, for R ^ Xj, the integral in ([70|) is 
roughly proportional to the string area RXj because Xj 
is an estimate for the string width: 



The e nerg ies ( [70| ) or (71) should be added to the esti- 
mate ( [43|) for the energy cost of the subsurface vortex 
termination. 

Taking, for example, Bi-2212 with s = 15 A, A,/ = 
4500 A, K = Xab/iab = 50, we have for i? = 1 /i: 



Ae 



167r2A 



In 



\K.Xai 



R 



(47r)3/2 



2R 

ttXj 



(72) 



We estimate the logarithm here as « — 2 (disregarding 
the factor (47r)~3/2 which could reduce this estimate to 
—4). The Josephson correction is 2R/ttXj w 1.4. Thus, 
in micron-size samples of Bi-2212, the termination of vor- 
tices under the surface is possible. 

The above treatment is concerned with the case of zero 
applied field. The energy cost of the subsurface termi- 
nation should increase with increasing applied field. To 
have an estimate of this increase, one can evaluate the 
work needed to remove a pancake from the top layer of a 
cylindrical samle of radius R in the presence of a small ex- 
ternal field Ha- Since one layer screening ability is weak, 
we assume the field uniform, which corresponds to the 
asymuthal current density j ~ cA/AttX'^ — cHar/8nX'^ 
(the vector potential can be taken a,s A — Har/2). The 
work of the Lorentz force to remove a pancake from the 
center is now readily obtained as « Ha4>oR^ /16ttA. To 
compare with the estimate (|72| ) we write the extra cost 
due to the applied field as 



HaTtR^ 



167r2A 



(73) 



Since the expression in parentheses of Eq. (72) is on the 
order one, we can say that the energy gain of the sub- 
surface termination is lost in fields Ha ^ 4'o/t^R'^ which 



gives Ha ^ 10 G for i? = 1 /x. This sets an approximate 
upper bound on fields in which the subsurface termina- 
tion may occur. 

We would like to mention in conclusion that, since con- 
figurations considered here and similar ones are not for- 
bidden by topology, they should be taken into account 
while studying fluctuations in layered weakly coupled 
materials. They are certainly of interest for physics of 
layered organic superconductors in which Ac may reach 
a tens-of- micron size at low temperatures (see, e.g., Ref. 
y) as well as for superconducting multilayers where the 
Josephson coupling can be reduced by varying the inter- 
layer spacing. 
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